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[57] ABSTRACT 

The described and improved multi-arm invention of this 
application presents three strategies for adaptive con- 
trol of cooperative multi-arm robots which coordinate 
control over a common load. In the position-position 
control strategy, the adaptive controllers ensure that 
the end-effector positions of both arms track desired 
trajectories in Cartesian space despite unknown time- 
varying interaction forces exerted through a load. In 
the position-hybrid control strategy, the adaptive con- 
troller of one arm controls end-effector motions in the 
free directions and applied forces in the constraint di- 
rections; while the adaptive controller of the other arm 
ensures that the end-effector tracks desired position 
trajectories. In the hybrid-hybrid control strategy, the 
adaptive controllers ensure that both end-effectors 
track reference position trajectories while simulta- 
neously applying desired forces on the load. In all three 
control strategies, the cross-coupling effects between 
the arms are treated as “disturbances” which are com- 
pensated for by the adaptive controllers while follow- 
ing desired commands in a common frame of reference. 
The adaptive controllers do not require the complex 
mathematical model of the arm dynamics or any knowl- 
edge of the arm dynamic parameters or the load param- 
eters such as mass and stiffness. Circuits in the adaptive 
feedback and feedforward controllers are varied by 
novel adaptation laws. 
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METHOD AND APPARATUS FOR ADAPTIVE 
FORCE AND POSITION CONTROL OF 
MANIPULATORS 

5 

ORIGIN OF THE INVENTION 

The invention described herein was made in the per- 
formance of work under a NASA contract and is sub- 
ject to the provisions of Public Law 96517 (35 USC 202) 
in which the contractor has elected to retain title. 

This is a continuation of application Ser. No. 
07/638,580, filed on Jan. 7, 1991, now abandoned, 
which is a division of application Ser. No. 07/253,510, 
filed Sep. 30, 1988, which is a continuation-in-part of a ^ 
application assigned to the same assignee and identified 
as a METHOD AND APPARATUS FOR ADAPT- 
IVE FORCE AND POSITION CONTROL OF MA- 
NIPULATORS, filed on Apr. 6, 1987 having Ser. No. 
07/035,061, now U.S. Pat No. 4,860,215. 2Q 

BACKGROUND OF THE INVENTION 

1. Field of the Invention 

This invention relates to control systems for robotic 
manipulators and more particularly relates to adaptive 2 5 
control for a plurality of robotic arms, e.g. a multi-arm 
system in which the cross-coupling between robotic 
arms through a common load is treated as though it 
were an externally caused disturbance. 

2. Description of the Prior Art 30 

An adaptive control system invented by this inventor 

is described and claimed in the above-identified applica- 
tion and a considerable number of publications are set 
forth therein, which publications are incorporated 
herein by reference as though they were set forth ex- 35 
pressly herein. The prior art and references cited 
therein are called to the attention of the Patent Office as 
being of background relevance to this invention. As of 
this filing date, the identified application has received an 
Office Action and the following patents have been 40 
cited; Horack U.S. Pat. No. 4,547,858, Koyama et al 
U.S. Pat. No. 4,580,229; Sugimoto et al U.S. Pat. No. 
4,621,332; and Osuka U.S. Pat. No. 4,725,942. These 
patents do not have any significant relevance to the 
invention herein described and claimed. 45 

The system of the identified application is improved 
upon by being extended from a single arm approach to 
a dual-arm (or higher) approach with a novel method 
and apparatus being described for overcoming the 
cross-coupling that exists in a common load being ma- 50 
nipulated by independently controlled manipulators, or 
arms, in a multi-arm system. Additional background 
material relevant to the development and a fuller under- 
standing of this invention is given in the following para- 
graphs. 55 

During the past decade, robot manipulators (“arms”) 
have been utilized in industry for performing simple 
tasks, and it is foreseen that in the near future anthropo- 
morphic robots will replace human operators in carry- 
ing out various complex tasks both in industry and in 60 
hazardous environments. Nevertheless, present-day 
robots can be considered at best as “handicapped” oper- 
ators due to their single-arm structure. It is evident that 
a multiplicity of robot arms yields greater dexterity and 
increased efficiency and provides capability of handling 65 
larger loads. Dual-arm robots will therefore have capa- 
bilities which may match those of ambidextrous human 
operators in dexterity and efficiency. 
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The research on dual-arm robots is at its early stages 
at the present time and a few approaches are currently 
available. E. Nakano et al., Cooperational Control of the 
Anthropomorphous Manipulator MEL ARM, Proc. 4th 
Intern. Conf. on Industrial Robots, pp. 251-260, 1974, 
propose a method for control of dual-arm robots in a 
master/slave manner. T. Ishida, Force Control in Coordi- 
nation of Two Arms , Proc. 5th Intern. Conf. on Artificial 
Intelligence, pp. 717-722, 1977, considers parallel and 
rotational transfer of loads using dual-arm robots. S. 
Fujii et al.. Coordinated Computer Control of a Pair of 
Manipulators, Proc. 4th World Congress on Theory of 
Machines and Mechanisms, pp. 411-417, Newcastle- 
upon-Tyne, England, 1975, suggest a technique for 
dual-arm control based on the method of virtual refer- 
ence. C. Alford et al., Coordinated Control of Two Robot 
Arms, Proc. Intern. Conf. on Robotics, pp. 468-473, 
Atlanta, Ga., 1984, describe a method for coordinated 
control of two arms. Y. Zheng et al., Constrained Rela- 
tions Between Two Coordinated Industrial Robots , Proc. 
Machine Intelligence Conf., Rochester, N.Y., 1985 and 
Computation of Input Generalized Forces for Robots with 
Closed Kinematic Chain Mechanisms , IEEE Journal of 
Robotics and Automation, pp. 95-103, Vol. RA-1, No. 
2, 1985, obtain constrained relations and control laws 
for two coordinated arms. T. Tam et al., Coordinated 
Control of Two Robot Arms, Proc. IEEE Intern. Conf. 
on Robotics and Automation, pp. 1193-1202, San Fran- 
cisco, Calif., 1986, employ the “Global” linearization 
technique for dual-arm control. S. Hayati, Hybrid Posi- 
tion/Force Control of Multi-Arm Cooperating Robots, 
Proc. IEEE Intern. Conf. on Robotics and Automation, 
pp. 82-89, San Francisco, Calif., 1986, and in pending 
U.S. Patent Application Method and Apparatus for Hy- 
brid Position /Control of Multi-Arm Cooperating Robot, 
filed Mar. 21, 1988, Ser. No. 06/845,991 proposes a 
method for controlling dual-arm robots based on parti- 
tioning the load between the arms. A. Koivo, Adaptive 
Position-Velocity-Force Control of Two Manipulators, 
Proc. 24th IEEE Conf. on Decision and Control, pp. 
1529-1532, Ft. Lauderdale, Fla., 1985, suggests an 
adaptive control technique for dual-arm robots using 
the self-tuning approach. J. Lim et al., On a Control 
Scheme for Two Cooperating Robot Arms, Proc. 24th 
IEEE Conf. on Decision and Control, pp. 334-337, Ft. 
Lauderdale, Fla., 1985, describes a positional control 
scheme for two cooperating robot arms. 

Some recent results of single-arm adaptive control 
are reported in the following papers: H. Seraji, Adaptive 
Control of Robotic Manipulators, JPL Engineering 
Memorandum 347-182, January, 1986; H. Seraji, Direct 
Adaptive Control of Manipulators in Cartesian Space, 
Journal of Robotic Systems, February, 1987 (to appear); 
and H. Seraji, Adaptive Forces and Position Control of 
Manipulators, JPL Engineering Memorandum 347-192, 
October, 1986. 

The above-identified articles, to the extent that they 
are properly considered prior art, do not teach or sug- 
gest that a dual-arm adaptive control system, nor such a 
system having adaptive hybrid control of each arm 
independently. Moreover, only in this application is it 
taught that a multi-arm adaptive control system is reli- 
ably operable so long as the load’s inter-arm cross-cou- 
pling is treated as though that cross-coupling were an 
externally caused disturbance. The adaptive hybrid 
control system of this invention can compensate for that 
cross-coupling because of the novel force and/or posi- 
tion control laws as herein defined. 
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There are certain key differences between my single- 
arm invention’s approach and the conventional hybrid 
control approach of Raibert and Craig, referred to 
above. Firstly, in my single-arm invention, the force or 
position control problems are formulated in the Carte- 
sian space with the end-effector Cartesian forces as the 
manipulated variables; whereas in Raibert and Craig, 
the problems are formulated in the joint space. The 
single-arm invention’s formulation results in computa- 
tional improvement since inverse Jacobians are not 
required for the controllers’ operation. Secondly, the 
single-arm invention’s hybrid system operates on the 
measured variables so as to produce the position and 
force variables that need to be controlled; whereas in 
Raibert and Craig, a selection matrix and its comple- 
ment are used after formulation of tracking errors. In 
summary, the Raibert and Craig disclosure simply does 
not teach or suggest the novel concepts of my inven- 
tion. 

SUMMARY OF THE INVENTION 

My invention of the parent application (“single-arm” 
invention) discloses a novel method and apparatus for 
the design of adaptive force and position controllers 
within a hybrid control architecture. The hybrid con- 
troller includes an adaptive force controller which 
achieves tracking of desired force setpoints, and an 
adaptive position controller which accomplishes track- 
ing of desired position trajectories. Force and position 
controller gains are varied by novel adaptation laws. 
These controllers are computationally fast and suitable 
for on-line implementation with high sampling rates and 
such adaptive gains are implemented in feedback and 
feedforward controllers. These controllers are capable 
of compensating for dynamic cross-couplings that exist 
between position and force control loops. 

A key feature of my single-arm invention is that my 
adaptive hybrid controller architecture does not require 
knowledge of the complex dynamic model or parameter 
values of the manipulator or the environment. The 
force and position controllers are linear and stable and 
generate real-time signals which vary and compensate 
for system non-linearities in order to achieve a desired 
position/force response. 

An additionally attractive feature of my single-arm 
adaptive controllers is an ability to compensate for dy- 
namic cross-couplings that exist between the position 
and force control loops in the hybrid control architec- 
ture. Furthermore, the adaptive force and position con- 
trollers have “learning capabilities” to cope with unpre- 
dictable changes in the manipulator or environment 
parameters such as environment or robotic arm stiff- 
ness. This is due to the fact that the controller gains are 
adapted rapidly on the basis of the manipulator’s Carte- 
sian space performance. Low computational require- 
ments make the control loops of the single-arm inven- 
tion suitable for implementation in on-line hybrid con- 
trol with high sampling rates. 

The described and improved multi-arm invention of 
this application presents three strategies for adaptive 
control of cooperative dual-arm robots. In the position- 
position control strategy, the adaptive controllers en- 
sure that the end-effector positions of both arms track 
desired trajectories in Cartesian space despite unknown 
time-varying interaction forces exerted through the 
load. In the position-hybrid control strategy, the adapt- 
ive controller of one arm controls end-effector motions 
in the free directions and applied forces in the constraint 
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directions; while the adaptive controller of the other 
arm ensures that the end-effector tracks desired position 
trajectories. In the hybrid-hybrid control strategy, the 
adaptive controllers ensure that both end-effectors 
track reference position trajectories while simulta- 
neously applying desired forces on the load. 

In all three control strategies, the cross-coupling 
effects between the arms are treated as “disturbances” 
which are rejected, or compensated for, by the adaptive 
controllers while following desired commands in a 
common frame of reference. The adaptive controllers 
do not require the complex mathematical model of the 
arm dynamics or any knowledge of the arm dynamic 
parameters or the load parameters such as mass and 
stiffness. The controllers have simple structures and are 
computationally fast for on-line implementation with 
high sampling rates. Simulation results are given to 
illustrate the proposed adaptive control strategies. 

BRIEF DESCRIPTION OF THE DRAWING 

FIG. 1 depicts a block diagram of an improved hy- 
brid control architecture in keeping with this invention; 

FIG. 2 is a simplified block diagram of a controller 
for performing adaptive force control in accordance 
with the invention; 

FIG. 3 is a simplified block diagram of a controller 
for performing adaptive position control in accordance 
with the invention; 

FIG. 4 is a block diagram of a basic module which is 
repeated as necessary for implementation of the control- 
ler functions of adaptive force and position control in 
accordance with this invention; 

FIG. 5 depicts a tri-level hierarchal control for two, 
or higher-numbered arm, robotic systems; 

FIGS. 6 and 7 depict, respectively a single-arm and a 
cooperative dual-arm manipulator; 

FIG. 8 depicts a simplified load diagram that is useful 
in developing a better understanding of the invention; 

FIGS. 9a, 9b, 9c, 9 d, 9e and 9 f depict the results of a 
computer simulation in adaptive position tracking 
which are useful in performance evaluation of the in- 
vention; 

FIGS. 10 a, 10 b, 10c, IQd and 10c, depict the results of 
a computer simulation in adaptive position-hybrid con- 
trol which are useful in performance evaluation of the 
invention; and 

FIGS. 11a, 116, 11c and 11 d, depict the results of a 
computer simulation in adaptive hybrid-hybrid control 
which are useful in performance evaluation of the in- 
vention. 

DESCRIPTION OF THE PREFERRED 
EMBODIMENT 

The description of the invention claimed in this appli- 
cation will first include a summary of the novel single- 
arm force and position control laws and then those laws 
will be developed and applied to a multi-arm system of 
this invention. Incorporated herein by reference, in 
accordance with § 608.01 (p) of the Manual of Patent 
Examining Procedure, from the above-identified single- 
arm application, is the mathematical formulation of the 
adaptive force and/or position control laws as set forth 
in section 1.1, page 10 through section 3, page 28. More 
particularly such incorporation includes Equations (21) 
through (24) for the gains of the PID controller of FIG. 
3 hereof in the implementation of the linear adaptive 
force control law given by Equations (25) and (26): and 
likewise the formulation and implementation of the 
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linear adaptive position control law of Equation (37) 
and the controller gains of Equations (38) through (43). 

To avoid confusion, the Equations of this continua- 
tion-in-part application are numbered starting with 
Equation no. 101, etc. The exact same force and posi- 5 
tion control laws apply to multi-arm hybrid adaptive 
control as apply to single-arm hybrid adaptive control 
except the mathematical notation is changed slightly. It 
will be readily apparent to those of ordinary skill in this 
art, however, that the Equations are the same and one 10 
of the keys to the inventive features of this application 
is noting that the cross-coupling terms, C p> Equation 
(25) and C/ Equation (37) in the adaptive position and 
force loops of a single-arm system, also compensates for 
the cross-coupling that the independently controlled 15 
arm of one controller introduces into another indepen- 
dently controlled arm in a dual or higher-numbered 
multi-arm system. Stated in a slightly different way the 
cross-coupling term which results from an environment 
in the single-arm disclosure is replaced by the coupling 20 
through the load that is being held by a second arm in 
the multi-arm disclosure. Otherwise, the two systems 
are the same and the analysis for one system is fully and 
completely applicable to the other system. 

The blanket generality noted above may be appreci- 25 
ated by a brief reference to FIG. 6 in comparison to 
FIG. 7. FIG. 6 depicts a common way of showing a 
robot manipulator 600 having at least two joint-con- 
nected sections 600A, 6O0B and an end effector 605. 
The manipulator 600’s end-effector 605 is mqveable in a 30 
Cartesian space 630 and may come into contact with 
another object or a wall, such as environment 610. In 
FIG. 7, two arms 600 and 660 are shown and the envi- 
ronment in this instance is replaced by common load 
610' that is being manipulated by and between the two 35 
arms, 600 and 660, respectively. 

In order to show the correspondence with my earlier 
invention, I have identified the common load by num- 
ber 610' so that its correspondence to the wall, or other 
object 610 will be more readily apparent. Although not 40 
depicted, common reference systems, independent con- 
trollers with synchronization at the trajectory generator 
level and common coordination (at the INTERMEDI- 
ATE LEVEL, FIG. 5) is to be understood throughout 
the disclosure as will become readily apparent by the 45 
following detailed description. 

The presentation of the detailed description given 
hereinafter in this application is structured as follows. In 
Section 1 (“position-position”), the position-position 
control strategy is discussed and documented by Simula- 50 
tion results. In Section 2 (“position-hybrid”), the posi- 
tion-hybrid, eg. position-force control strategy is devel- 
oped and demonstrated by a numerical example. The 
hybrid-hybrid, force/position, or vice-versa, control 
strategy is addressed in Section 3 (“hybrid-hybrid”) and 55 
an illustrative example is given. The remainder of this 
application discusses the results of the foregoing multi- 
arm presentation and draws some conclusions. 

SECTION 1— POSITION-POSITION 60 

In this section, we shall investigate the first control 
strategy for dual-arm manipulators in which both arms 
are in pure position control, as shown in FIG. 2. In 
other words, the positions and orientations of both end- 
effectors are required to track desired trajectories in a 65 
common frame of reference. In this situation, uncon- 
trolled forces and torques will be exerted on the com- 
mon load 610' held by the end-effectors 605,665 of the 
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dual arms such as those represented in the LOW 
LEVEL section of FIG. 5. Investigated herein is the 
performance of the adaptive position control law of my 
system in the face of the interaction forces and torques 
exerted through the load 610'. 

The dynamic model of each manipulator arm (arm 
600 or 660, FIG. 7) can be represented by a differential 
equation in Cartesian space as was observed by O. Kha- 
tib, in Dynamic Control of Manipulators in Cartesian 
Space , Proc. 6th IFToMM Congress on Theory of Ma- 
chines and Mechanisms, pp. 1128-1131, New Delhi, 
India, 1983. 

M{X)X+N(X,V+G{X)+H{X)±f=F ( 101 ) 

where the above terms are defined as: 

X,X,X=nX 1 vectors of end-effector position, veloc- 
ity and acceleration in a fixed task-related Carte- 
sian frame of reference 

F=nX 1 vector of “virtual” Cartesian forces applied 
to the end-effector as the control input 

M (X)=nXn symmetric positive-definite Cartesian 
mass matrix 

N(X,X)=nXl Cartesian Coriolis and centrifugal 
force vector 

G(X)=nX 1 Cartesian gravity loading vector 

H(X)=nXl Cartesian friction force vector 

f=NX 1 vector of forces and torques exerted by the 
end-effector on the load. 

In the above-noted example the load is the shared 
common load 610' between the dual-arm system of 
FIGS. 5 and 7. 

The force/torque vector f both imparts motion to and 
applies force/torque on the load 610' and acts as the 
coupling element between the two arms 600 and 660. In 
the following analysis, the force/torque vector f will be 
considered as a “disturbance input” (analogous to my 
earlier-described auxiliary input signal d(t), C p and C/) 
to the position control system. The function of the con- 
trol system is to ensure that the end-effector position 
vector X tracks the nXl vector of desired trajectory 
Xd despite the disturbance force f. For each manipulator 
arm, let us apply the linear adaptive position control 
law described in my earlier application or it’s corre- 
sponding paper which is identified as H. Seraji, Direct 
Adaptive Control of Manipulators in Cartesian Space , 
Journal of Robotic Systems, February, 1987. Such a 
position control law is: 

XM+AWU 01 ( 102 ) 

as shown in FIG. 3, where E(t)=X<f(t)— X(t) is the nX 1 
position tracking-error vector. In the control law (102), 
the nX 1 vector d(t) is an auxiliary signal to be synthe- 
sized by the adaption scheme, while [K^E+KyE] and 
[CX^+BX^+AX^] are the contributions due to the 
feedback and feedforward controllers respectively. Fol- 
lowing my earlier-described method, the required auxil- 
iary signal and controller gains are updated according 
to the following adaptation laws: 

r t ( 103 ) 

d(t) = d{ 0) + 8i Ht)dt + 6 2 *0 

J 0 

r t ( 104 ) 

Kp(t) = 0) + ai Ht)E{t)dt + a 

J 0 
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•continued 


K^i) = KM + 0i 

f ? KtWm + faKtWit) 
J 0 

(105) 

£ 

II 

O 

o 

4 

t 

Kwm + V2 «« 

0 

(106) 

B(t) = 5(0) + 71 f 

t 

rif)XJ{t)dt 4 72<W(0 

0 

(107) 

A(t) = A( 0) + \i J 

r KWP 4 *2*W(0 
0 

(108) 

where 

rit) = WpE(i) + Wfi(t) 

(109) 


in an nXl vector, {Si, a\, fi\, vi, 71, Ai} are positive 
scalars, {82, a2, £2, V2, 72. A2} are positive or zero sca- 
lars, and the prime denotes transposition. In equation 
(109), Wp and W v are nXn constant weighing matrices 
chosen by the designer to reflect the relative signifi- 
cance of the position and velocity errors E and E. It 
must be noted that since we cannot physically apply the 
Cartesian control force F to the end-effector, we instead 
compute the nX 1 equivalent joint torque vector T to 
effectively cause this force. Thus, for each manipulator 
arm (arm 600 or 660, FIG. 7, for example), the control 
law in joint space is given by 

7 X0 = -+-^(0^(0 -t- h- - 
C{t)X^B{f)X^A{t)X^t)} (1 10) 

where 6 is the nX 1 vector of joint angular positions and 
J(0) is the n X n Jacobian matrix of the manipulator arm. 

Because of the simplicity of the adaptation laws (103) 
through (108), the robot control algorithm can be im- 
plemented using high sampling rates (typically 1 KHz). 
In each sampling period (—1 msec), the controller gains 
can change significantly; whereas the terms M, N, G, H, 
and f in the robot model (101) cannot change notice- 
ably. As a result, in deriving equations (103) through 
(108), it was assumed that these terms are unknown and 
“slowly time-varying” relative to the gain variations as 
such are changed by these adaptation laws. It is seen 
that the inclusion of the disturbance force f in the robot 
model (101) does not affect the controller adaptation 
laws since the change in f over one sampling period is 
relatively small. 

The above-noted observation likewise suggests that 
when both manipulator arms are controlled using the 
two independent adaptive position controllers as imple- 
mented by my adaptation laws, we observe that the 
end-effectors will track the desired position trajectories 
despite the interaction forces and torques exerted 
through the load. It must be noted that since the force 
on the load is not controlled variable in this invention, 
this strategy can lead to undesirable load forces when 
the position trajectories are not planned in coordination 
or are not tracked closely. The position-position control 
strategy is illustrated by the following example. 

Consider the planar dual-arm manipulator in a hori- 
zontal plane (g=0) shown in FIG. 7, where the two 
arms 600 and 660 are identical and each arm has two 
links and an end-effector with negligible inertia. Sup- 
pose that the load is a linear spring of natural length \ 0 
and coefficient of stiffness (elasticity) K^. In the joint 


8 

space, the dynamic equations of motion for each arm 
which relates the joint torque vector 



10 to the joint angle vector 


0 = 


15 


is given by H. Seraji, et al., Linear Multivariable Control 
of Two-Link Robots , Journal of Robotic Systems, pp. 
349-365, Vol. 3, No. 4, 1986 as follows 

T(t)=M(0)0+N(0, 0)+H(0)±J'(0)f (1 1 1) 


where the above terms are defined as: 

25 Inertia matrix 


30 


35 


40 


45 


M(6) » 


a\ + G2COS02 a3(<22/2)cos02 ^ 
03 4 O2/2)cos02 a 3 


Coriolis and centrifugal torque vector 

( — (<32sin02) (0102 4 0l/2) 


N\ 0,0 = 


^(a2sin02)0i 2 /2 
Coulomb and viscous friction torque vector 


H{Q) * 


f V\b\ + K2Sgn(0O ^ 
^302 4- K4Sgn(0 2 ) ^ 


Jacobian matrix 
/ 


m = 


— /jsin0i — /2sin(0i -I- 02) — /2sin(0i + 02) 

/icos0i -f /2COS(0l 4 02) /2cos(0i + 02) 


In the above expressions, at, . . . , as are constant param- 
50 eters obtained from the masses {mi, m2} and the lengths 
{U, I2} of the links, {vi, V3} and {v2, V4} are coefficients 
of viscous and Coulomb frictions respectively, and the 
links are uniform with centers-of-gravity located at 
mid-lengths. The interaction force vector f is modelled 
55 as 


1 


( K 0 (l - l 0 )cosa > 

s 

II 

4/ 

1 K 0 (l — / 0 )sina 


where 1 = V (x2 — x 1) 2 + (y 2 “ y l) 2 is the instantaneous 
length of the coordinates of each end-effector given by 


x=do- f/i cos 0i 4 / 2 ( 0 1+02) 
y—l\ sin 0i 4/2 sin (01 + 02) 
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and do is the base distance from the origin and a is the 
angle between the spring and the x-axis, that is 


the terms in the control law (112) are adapted as fol- 
lows: 


yi-y\ 

tana = . 5 

x 2 - *1 J 

The mass of the spring is assumed to be negligible com- 
pared to the masses of the links and hence the force in 
the spring dominates over the inertial force required to 
move the spring. 

For the particular system under study, the numerical 
values of the parameters are chosen as 

mi =15.91 kg; m 2 = 11-36 kg; h=l 2 =-0.432 .m ^ 

V 1 = 6.0Nt.m/rad.sec“ \ V3 =4.0Nt.m/rad.sec - 1 

V 2 = V 4 = l.ONt.m; Ko=1000.0Nt-/m 

The natural length of the spring is set equal t o the initial 20 
distance bet ween the end-effectors; i.e. loV(x 2 /— x//) 2 " 
+(y2z— y//) 2 so that there will be no initial force in the 
spring. The masses and lengths of the links are those of 
links 2 and 3 of the Unimation PUMA 560 arm. It must 
be emphasized that the mathematical model and the 25 
parameter values of the arms and the spring are used 
merely to simulate the system behavior and are not used 
to generate the control action. This invention does not 
depend upon knowledge of the complex mathematical 
model of the arm dynamics or any knowledge of the 30 
arm dynamic parameters or the load parameters such as 
mass or stiffness. As was noted hereinbefore, the lack of 
requirement of such knowledge must be acknowledged 
as being one of the key features that contributes to the 
simplified nature of this adaptive hybrid system inven- 35 
tion and its novel force and/or position control laws, or 
algorithms. 

Suppose that the end-effector of the right arm is re- 
quired to track the desired position trajectories 

x r d( t ) = 0. 864[ 1 + 3 exp ( — r/0.3) — 4 exp ( - t/Q.4]m. 

JVd(0= 0 *864[3 exp {—t/ 0.3)— 4 exp (— r/0.4]m. 

to move from the initial point {0,0,-0.864} to the final 45 
point {0.864,0.0} on the straight line 
Y«f(t)=Xrd(t)— 0.864. The end-effector of the left arm 
is required to track the desired position trajectories 

x/</f)=2.0+0.432[l 4.3 exp (— f/ 0.3)— 4 exp 50 

(-«/ 0.4)]m. 

yid( t) = 0.432 [ — 1 + 3 exp ( — t/Q. 3) — 4 exp 
(— f/0.4)]m. 

to move from the initial point {2.0,-0.864} to the final 
point {2.432,-0.432} on the straight line 55 
ya(t)=xa(t)— 2.864. In this physical set-up the bases of 
the two arms are 2.0 meters apart, and the desired posi- 
tion trajectories for the right and left arms were chosen 
such that appreciable forces were created in the spring. 

The two arms, in accordance with my invention, are 
controlled by identical and independent adaptive posi- 
tion control laws given by 

aw) 

BWM+AW<tf)} ( 112 ) 65 

where E(t)=X^t)— X(t) is the position tracking-error, 
X(t)=[x(t),y(t)]', and X<f( r )= [*<#), y<*(t)]'. For both arms, 


d(t) = 0.5*0 -h 0.5 f * 
J 0 

K p (f) = 2 J 

r* t 
0 

KM = 2 J 

. t 

*t)E(i)d 

0 

d 

II 

C KWP 
} 0 

B(t) = 0.5 

f Kmm 

X 

II 

p 

Ul 

f ? *t)XMdt 
' 0 

where 

*0 = 8000E(0 + 800i(0 


Note that the initial values of the controller terms are all 
chosen arbitrarily as zero. A simple trapezoidal rule is 
used to compute the integrals in the adaptation laws 
with dt= 1 msec. 

To evaluate the performance of the position-position 
control strategy, the nonlinear dynamic models of the 
arms (Equation 111) and the linear adaptive control 
laws (Equation 112) are simulated on a DEC- VAX 
11/750 computer with the sampling period of 1 msec. 
The results of this simulation are shown in FIGS. 9(g) 
through 9(f). FIG. 9(g ) shows the magnitude of the 
interaction force f m = V f x 2 4* fy 2 acting on the end-effec- 
tors through the spring. It is seen that f m changes from 
the initial value of zero to the final value of 373 New- 
tons in 2 seconds. FIGS. 9(b) through 9(e) show that the 
end-effector coordinates x(t) and y(t) of both arms track 
their corresponding reference trajectories x^(t) and 
yd(t) very closely despite the large interaction force f m . 
The paths followed by the end-effectors in the horizon- 
tal plane are shown in FIG. 9(f). Note that both end- 
effectors track the straight lines as desired. 

We conclude that the adaptive controller for each 
arm performs remarkably well in faithful position track- 
ing despite very large unknown and time- varying inter- 
action forces between the arms. 

SECTION 2— POSITION-HYBRID 

In this section, the position-hybrid control strategy 
for dual-arm manipulators will be studied in which the 
left arm is in pure position control and the right arm is 
in hybrid position/force control, as shown in one em- 
bodiment of FIG. 3. In other words, for the left arm 
600, FIG. 7, the end-effector position is required to 
track a desired trajectory in a frame of reference. For 
the right arm 660, FIG. 7, in the same reference frame, 
the contact force between it’s end-effector 665 and the 
load 610" must be controlled in the directions con- 
strained by the load 61 0 f , while the end-effector 665’ s 
position is to be controlled simultaneously in the free 
directions. This control strategy is also applicable when 
one robot arm is confined for operation only in a posi- 
tion control mode whereas, simultaneously, the other 
arm can be controlled in hybrid control mode. 
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For the left arm 600, the interaction forces and 
torques exerted through the load 610' are considered as 
“disturbances,” and the adaptive position control sys- 
tem can ensure tracking of the desired position trajecto- 
ries despite such disturbances, as outlined in Section 1 . 
The adaptive position control law for the left arm 
shown in FIG. 7 is given by Equation (37) in my earlier 
application and may be applied here as 

(113) 

where T/ is the nX 1 joint torque vector, 0/is the nX 1 
joint angle vector, J /(#/) is the nXn Jacobian matrix, 
X/d(t) — X(t) is the nXl position tracking-error vector 
and the terms in Equation (113) are adapted as follows: 


d(t) = 7(0) + 81 J 

t 

r{t)dt -f- 8 2 r(t) 

0 

(114) 

r/r) = Kp(0) + a. 

f r + a 2 r{t)E{t) 

J 0 

(115) 

m = ka o) + Hi 

f ? rit)E{t)dt + p 2 Ht)E{t) 
J 0 

(116) 

C(t) = C( 0 ) + VI J 

r r(t)XiJ(t)dt + V 2 r®XM 
0 

(117) 

Bit) = 1(0) + yi J 

t 

+ 72 
0 

(118) 

A{t) = A( 0) + \i J 

t 

r{i)X ld {t)dt + \2*t)Xldt) 

0 

(119) 

K0== WpEii) + W r E(t) 

(120) 


and the symbols are defined in Section 1. 

We shall now discuss the hybrid position/force con- 
troller for the right arm 660, FIG. 7. Consider a task- 
related “constraint frame” (coordinate system) which is 
defined by the particular contact situation occurring 
between the right end-effector 665 and the load 610'. In 
this frame, the n degrees-of-freedom (or directions) in 
the Cartesian space {X} can be partitioned into two 
orthogonal sets; the m constraint directions in subspace 
{Z} and the 1 free directions in subspace {Y}, with 
n=m+l. In the m constraint directions, the end-effec- 
tor makes contact with the load 610' and the contact 
force needs to be controlled. In the 1 free directions, the 
end-effector is free to move and the end-effector posi- 
tion is to be controlled. In the hybrid control architec- 
ture (as noted in the aforesaid Raibert and Craig and 
Mason articles), two separate controllers may be em- 
ployed for simultaneous force and position control. The 
“force controller” achieves tracking of desired force 
setpoints in the constraint directions; while the “posi- 
tion controller” accomplishes tracking of desired posi- 
tion trajectories in the free directions. 

The dynamic model of the right arm in the constraint 
directions can be written as 

A(X,X)P(t) +B(X,X)P(t) + F®+ Cp( Y)±f z =F&) ( 121 ) 

where F z is the mXl “virtual” Cartesian force vector 
applied to the end-effector in the constraint directions, 
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Z is the mX 1 vector of end-effector position, the mXm 
matrices A and B are highly complex nonlinear func- 
tions of the end-effector position X, C p is the cross-cou- 
pling from the position loop into the force loop and f z is 
5 the component of the force exerted on the end-effector 
by the load 610' in the constraint directions. The term f z 
represents the cross-coupling that exists between the 
arms 600,660 through the load 610' and is considered as 
a “disturbance” to the hybrid controller. 

10 An adaptive force control scheme is developed 
within the hybrid control architecture as claimed 
herein. For the right arm, 660, the linear adaptive force 
control law in the constraint directions is given by 

Equation 122 below as 
15 

+ d(t\+ K p {t)E{t) + Kj(t) f o- 

'EiOdt-KvCOZtf ( 122 ) 

where P j(f) is the desired contact force on the load 610' 
20 used as a feedforward term, d(t) is an auxiliary signal, 
E(t)=Pr(t)— P(t) is the deviation of the actual force P(t) 
from the desired value, and {Kp(t) , Kj(t) , Kv(t) } are 
adaptive gains of the PID controller. The terms in the 
force control law. Equation (122), are adapted as fol- 
25 lows: 



d( t ) = d( o) + 8 , r r 
j 0 

q(t)dt + S 2 q(t) 

(123) 

30 

KKt) = KK 0 ) + a, J 

. t 

4 WW + a 2 g(t) E*{t) 

0 

(124) 

35 

K p (t) = K p i 0) + Hi J 

! ‘ qmit)dt + thq(')m 
0 

(125) 


KAO = KAO) + 71 j 

• t 

q{t)Z'{t)dt + 72 q(t)Z{t) 

0 

(126) 

40 

where 

<7(0 - m E*(t) + WpE{t) - WyZ{t) 

(127) 


In equations (123) through (127), 

45 t 

B\t) = E{t)dt 
J 0 

is the integral error vector, {Si,ai,/?i,yi} are positive 
50 scalars, { 82 , 02 ,^ 2 , 72 } are positive or zero scalars, and 
{W/,Wp,W v } are constant weighting matrices chosen 
by the designer to reflect the relative significance of E*, 
E and Z. 

The dynamic model of the right arm in the free direc- 
55 tions can be written as 

AJiX,X)Y(t)+B 0 (X,X)Y(t)+ C 0 (X,X) Y(t) + Cj(P)- 

+fy=F>tf) (128) 

where f y is the component of the end-effector force in 
the free directions, C/is the cross-coupling from the 
force loop, A 0 , B& C 0 are complex nonlinear matrices, Y 
is the end-effector position vector and Fj,is the “virtual” 
end-effector control force. For the right arm, the linear 
adaptive position control law in the free directions is 
65 given by 

F/t)=^t)±K^t)E^K^ • .. 

( 129 ) 
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as in Section 2, where R is the desired position trajec- 
tory, Bp—R — Y is the position tracking-error, and F}is 
the “virtual” Cartesian force in the free directions. 
Thus, in order to implement the force and position 5 
controllers in accordance with Equations (122) and 
(129) in the hybrid control architecture, the joint space 
control law for the right arm is given by 



where 0 r is the joint angle vector, T r is the joint torque 15 
vector, and J r is the Jacobian matrix of the right arm 
with appropriate reordering of columns of J r if neces- 
sary. 

The hybrid controller adaptation laws, Equations 
(103) through (108) and (123) through (126) are ex- 20 
tremely simple, and therefore the control algorithm can 
be implemented using high sampling rates (rates (~1 
KHz); yielding improved performance particularly in 
force control applications. Since in each sampling per- 
iod (~lm sec) the terms in the robot models (121) and 25 
(128) cannot change noticeably, it is reasonable to as- 
sume that these terms are “slowly time-varying” com- 
pared to the adaptation scheme. Thus the inclusion of 
the disturbance f 2 and f^in the robot models (121) and 
(128) does not effect the controller performance. 30 

It is concluded, therefor, that using the position- 
hybrid control strategy, the left end-effector will track 
the desired position trajectory despite the interaction 
forces through the load. The right end-effector will 
exert the desired force on the load in certain directions 35 
while simultaneously tracking the desired position tra- 
jectory in the orthogonal directions. It must be noted 
that in this control strategy, slight fluctuations may be 
observed on the load force due to very small vibrations 
of the left arm under position control. The position- 40 
hybrid control strategy is now illustrated by the follow- 
ing example. 

Let us consider the dual-arm manipulator discussed in 
Example 2 (with negligible friction) and suppose that 
the end-effectors carry a rigid load in a horizontal plane, 45 
as shown in FIG. 7. Suppose that the load 610', as 
shown, is initially parallel to the x-axis and we wish to 
make a parallel translation of the load 610' while the 
right arm, 660, applies a constant desired contact force 
on the load. This problem is similar to pulling out a 50 
drawer using two arms. 

The x-component of the force on the load 610' can be 
written as 

fx=k 0 {{x2~xi)-l 0 ] 55 


Suppose that the end-effector 605 of the left arm 600 
is required to track the desired position trajectories 

0.432 meter 

TVXO— 0-432[l +6 exp (—r/0.3) — 8 exp 
(— r/0.4)]meter 

to move from the initial point {0.432, —0.432} to the 
final point {0.432, 0.432} on the vertical line 
x/c^t)— 0.432. The end-effector 665 of the right arm 660 
is required to apply a specified force setpoint P r ^ on the 
load 610' in the x direction and simultaneously track the 
desired position trajectory y r d(t) in the y-direction; 
where 

P r d= 10 Newton 

JWO)«0432[l+6 exp (-r/0.3)- 8 exp (-r/0.4)] 
meter 

so that the end-effector moves from the initial point 
{2.432, 0.432} to the final point {2.432, 0.432} on the 
vertical line x r d(r)— 2.432, where 2.0 meter is the dis- 
tance between the bases of the two arms. 

The adaptive position control law for the left arm is 
given by 

+B{t)Xid!)+A(t)XiM (131) 

where E(t)=X/^(t)— X^(t) is the position tracking-error, 
Xi(t)=[x/(t), yKt)]', Xa(t)«pCa(t), y w(t)]' and the terms 


in equation (131) are: 

d(t) = 0.Sr(t) + 0*5 f 7{t)dt 
J 0 

04 

1! 

is? 

xtwm 

0 


t 

nt)E(t)dt 

0 

m - o.5 j 

0 

m = o.5 j 

'*mtkMdt 

0 

AM = 0.5 J 

. t 

HWiMdt 

0 

where 

7(0 = 2000QE(f) + 2000£(r) 


where (x 2 — xi) is the distance between the end-effec- 
tors, \ 0 is the length of the load, K 0 is the stiffness of 
force/torque sensors mounted on the end-effectors to 
measure the force f*. The y-component of the force on 
the load is 


Note that the initial values of the controller terms are 
set to zero. The integrals in the control law (131) are 
computed using a simple trapezoidal rule with dt=0.5 
msec. 

The adaptive hybrid position/force control law for 
the right arm is given by 


60 


r m .. 

fy = Y y 

65 

where the arms are assumed to share the load equally; 
i.e. each arm “sees” the load mass as m/2. In the simula- 
tion, we take K^= 1000.0 Nt/m and m= 1.0 Kg. 


7X0 - ftfr) 



\ 




(132) 
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where 

FAO = PM + <K0 + m)e x *(t) + 

W = 5(0 + + KMkyHt) + 

C( 0 ^r) + S(t)yA0 + mm 

eAO = P^r) — P^r) — force tracking-error 
ey(t) = »^/) — = position tracking-error 

e x *(r) = e x t)dt = integral force error 
J 0 

The adaptation laws for the force controller are: 

d(t) = AO + C q{t)dt 
J 0 


KAO =100+100 


> f * 

J 0 

K p {t) = 100 + 10 f ? 9(0*^ 

J n 


A'XO = 5000 - 40000 


J o 


where 

9(0 = lOe^) - 100 x (0 + 10^*(0 
The adaptation laws for the position controller are: 


10 


15 


20 


25 


30 


35 


d { 0 = 0 . 5 <t) + 0.5 j AOdt 
J 0 


K p (t) = 2 J 

. t 

K 0 e/ 0 dt 

0 

40 

m = 2 j 

. t 

mAOdt 

0 


C(t ) = 0.5 

f KOyAOdt 
1 0 

45 

V> 

O 

II 

/ *5 

f ^ KOMOdt 
! 0 

50 

A(t) = 0.5 J 

f r Komdt 
f 0 


where 


55 

AO = lOOOe/?) + 500 eAO 



The integrals in the above adaptation laws are evalu- 
ated using the trapezoidal rule with dt=0.5 msec. 

To evaluate the performance of the position-hybrid 60 
control strategy, the nonlinear dynamic models of the 
arms (Equation 111) and the linear adaptive control 
laws (Equations 131 through 132) are simulated on a 
DEC- VAX 11/750 computer with the sampling period 
of 0.5 msec and the simulation results are shown in 65 
FIGS. 10(a) through 10(£). These Figures show that the 
end-effector coordinates xj(t) and y/(t) of the left arm 
track the desired position trajectories x/^t) and yidff) 
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very closely. FIGS. 10(c) through 10(cT) indicate that 
the right end-effector exerts the desired force of 10 
Newtons on the load in the x-direction and tracks the 
desired position trajectory y r ^t) in the y-direction. No- 
tice that the oscillations in the force response in FIG. 
10(c) are due to the very small variations of the left 
end-effector x-coordinate as shown in FIG. 10(a). 
These oscillations can be reduced by increasing the 
gains of the position controller for the left arm. The 
paths traversed by the end-effectors in the horizontal 
plane are shown in FIG. 10(c). It is seen that the left 
end-effector tracks a vertical straight line, whereas the 
right end-effector moves in initially to produce the 
desired contact force and then tracks a vertical straight 
line. Thus, the adaptive position and hybrid controllers 
for each arm perform well in the dual-arm situation. 

SECTION 3— HYBRID-HYBRID 

In this section, the hybrid-hybrid control strategy for 
dual-arm manipulators will be studied in which both 
arms are in hybrid position/force control. In other 
words, in a common frame of reference for both arms, 
the forces exerted by the end-effectors on the load in the 
constraint directions {Z} must be controlled; while 
simultaneously the end-effectors are required to track 
desired position trajectories in the free directions {Y}. 
Any unwanted forces and torques on the load generated 
by the relative position and orientation of the end-effec- 
tors will act as “disturbances” and the adaptive hybrid 
controllers ensure that the desired position/force trajec- 
tories are tracked despite such disturbances. 

Following Section 2, for each manipulator arm the 
hybrid position/force control law in the joint space can 
be written as 


no = m 


'fao 

FA0 


\ 


J 


( 133 ) 


where J (6) is the Jacobian matrix (with appropriate 
column reordering if necessary), and F^t) and F^t) are 
the “virtual” Cartesian forces applied to the end-effec- 
tor in the constraint directions {Z} and free directions 
{Y}, respectively. The force control law is given by 


r t 034) 

FA!) = m + d{t) + KAO EAOdt + e^OEAO - KAOAO 
J 0 

where P^t) is the desired force setpoint, 
E^t) = P^t) — Pz(t) is the force tracking-error and the 
adaptation laws are: 


d{t) = d{ 0) + Si f f q{0dt + 5 2 g(0 
J o 

KAO - KAO) + ai f r <A0Ez*X0dt + a 2 <A0Ez*'(0 
J 0 

KAO = K/ 0 ) + 0 i f f q( 0 E z '( 0 dt + PlA 0 E z ( (0 
J 0 

= KAO) -yiC <A0Z{0dt - yiqiom 
J n 


KAO 
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-continued 

where 

*0 = W X E z *{t) + WpE^t) - WyZ(t ) and 

Ez*(t) = f ( £,(0* and {>£/, ^ W v } 

J 0 


are desired weighting matrices. 

The position control law is expressed as 

F/0=^(0+^>(OE/t)+EA0E/0+C(t)+mE(!- 

)+mm 
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{[0.432, 0.432], [2.432, 0.432]} while exerting a force of 
10 NT on the load. 

The two arms are controlled by identical and inde- 
pendent adaptive hybrid control laws given by 


10 


m = SO) 


''em ') 

) 


where 


(136) 


(135) F^t) = P r &) + d{t) + Kj{t)e x *{t) + K p {t)e^t) - K^t)x(t) 


where R(t) is the desired position trajectory, 
E/t)=R(t)— Y(T) is the position tracking-error and the 
adaptation laws are: 


7(0 = d( 0) + Si Ht)dt + fi 2 *0 
J 0 


K p {t) = K p { 0) + v\ iif)Ey(t)dt + i m)E y \t) 
J 0 


F^t) = 7(0 + A^O&fO) 4- K^t) + 

c(omo + my &) 

e^t) = P r d(0 — Pj&O — force tracking error 
e/0 — jPr/O — yAO = position tracking-error 

^*(r) = e/0^ = integral force error 

J 0 


«0 = *v(0) + Ti f ? Mfyw + m*t)EyV) 
J o 

C(r) - C(0) + jTi[ r r{t)R\t)dt + wm'to 
J n 


Bit) = 5(0) + 71 f ? KOH'Cd* + yiKt)R\t) 
J n 


A{t) = +(0) -h Xi f rit)R\i)dt + X 2 *0#(0 
J 0 

where 


*0 = WpEyAt) + ^i/0 


25 


30 


35 


40 


and {Wp, W v } are desired weighting matrices. 

The above controller adaptation laws are extremely 
simple and therefore the hybrid control algorithm can 
be implemented using high sampling rates (—1 KHz); 
yielding improved performance. Under the adaptive 45 
hybrid controllers, both end-effectors are expected to 
exert the desired forces on the load while simulta- 
neously moving on desired trajectories. The hybrid- 
hybrid control strategy is most available when simulta- 
neous control of both position and force is required. 50 
The following example illustrates the hybrid-hybrid 
control strategy. 

Consider the dual-arm manipulator and load, as dis- 
cussed above in conjunction with FIG. 7, and suppose 
that both end-effectors are required to exert a constant 55 
desired force on the load in the x-direction while mov- 
ing the load in the y-direction. In this situation, the x 
and y components of the load force are the same as 
developed in Section 2. 

Suppose further that the desired force and position 60 
trajectories for both arms are specified as * 


P jc/0 — ± 1 (5/Vewton 

y/0=0. 432[1 +6 exp (-r/0.3)-8 exp i~ t/0A)] meter 


65 


The adaptation laws for the force controller are: 

d(t) = q{t) + C qit)dt 
J 0 

Kjit) = 100 + 10 f 1 q{t)e x *{i)dt 
j 0 

K p it) = 100 + 4 f ? qQeMdt 
J 0 

KAt) = 5000 - 1000 f qit)xit)dt 
J 0 

where 

«0) = e*<0 - 10*0 + e x *(t) 

The adaptation laws for the position controller are: 

7(0 = 0.5r(0 +0.5 f r Kt)dt 
J 0 

K^t) = 2 f ? lifieffldt 

J 0 

KAt) = 2 C r(t)ey(t)dt 
J 0 

C(0 = 0.5 f 1 r{t)yAt)dt 
J 0 

Bit) = 0.5 f r rit)yAt)dt 
J 0 

Ait) — 0.5 f r r{t)yAt)dt 
J 0 

where 


so that the end-effectors move from the initial points 
{[0.432, -0.432], [2.432, -0.432]} to the final points 


K0 - SOOOe/t) + 2000e/0 
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The integrals in the above adaptation laws are evaluated 
using the trapezoidal rule with dt=0.5 msec. 

To evaluate the performance of the hybrid-hybrid 
control strategy, the nonlinear dynamic models of the 
arms (Equation 111) and the linear adaptive control 5 
laws (Equation 136) are simulated on a DEC- VAX 
1 1/750 computer with the sampling period of 0.5 msec 
and the simulation results are shown in FIGS. ll(a)-{d). 
FIG. 11(c) shows that the force exerted by the end- 
effectors on the load in the x-direction is equal to the 10 
desired setpoint of 10 Nt. FIGS. 11(6) through 11(c) 
indicate that the coordinates of the end-effectors in the 
y-direction track the desired position trajectories. The 
paths traced by the end-effectors in the horizontal plane 
are shown in FIG. 11(d)- It is seen that both end-effec- 15 
tors move in initially to produce the desired force and 
then track vertical straight lines. Thus the adaptive 
hybrid controllers perform extremely well for simulta- 
neous force and position control. 

Three adaptive control strategies for cooperative 20 
dual-arm robots have been described. In these strate- 
gies, each robot arm is considered a subsystem of the 
total system and is controlled independently using an 
adaptive controller in the low level of the control hier- 
archy. Each controller ensures that the controlled van- 25 
ables follow desired commands and reject unwanted 
cross-coupling effects from other subsystems which are 
treated as “disturbances.” The subsystems are coordi- 
nated through trajectory generators in the intermediate 
level, where synchronous desired trajectories for both 30 
arms are specified in a common task-related frame of 
reference. An important feature of the present approach 
is that the overall control system for N cooperative 
arms is reduced to N decentralized independent single- 
arm controllers. The control schemes do not require 35 
communication and data exchange among individual 
controllers, which is an appealing feature from both 
computational and reliability points of view. Further- 
more, available techniques for single-arm control can be 
utilized directly in multiple-arm environments. 40 

The control strategies described herein do not require 
the knowledge of the load parameters such as mass and 
stiffness or the robot dynamic parameters such as link 
masses and inertias, and can therefore cope with uncer- 
tainties or variations in the system parameters. Further- 45 
more, the complex dynamic model of the arms are not 
used in generating the control actions. The control 
schemes are very simple and extremely fast for on-line 
implementation with high sampling rates, yielding im- 
proved dynamic performance. The control methodol- 50 
ogy described herein can also be utilized in the coordi- 
nated control of N-arm robots when N exceeds two. 

APPENDIX 

In this Appendix, we obtain a simple expression for 55 
the force on a rigid load held by two end-effectors in a 
horizontal plane. 

Let us consider a rigid load of length loheld firmly by 
two end-effectors equipped with force/torque sensors 
as shown in FIGS. 7 and 8. The sensors are modelled as 60 
linear springs with stiffness coefficients Ki and K 2 and 
natural lengths li and I 2 . The forces Fi and F 2 exerted by 
the sensors on the rigid load are given by 

Fi=Ki(l'i—!i) (137) 65 
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where l'i and l^are instantaneous lengths of the springs. 
At equilibrium, the forces Fi and F 2 must be equal; 
otherwise the load will reposition itself under the net 
force to reach the equilibrium condition. Let us denote 
the force exerted by the springs on the load by 
F==Fi=F 2 ; hence 

F=F 1 =^ 1 [/' I -/ 1 -f/ 0 +/o + /' 2 -/' 2 +/ 2 -/ 2 ] (139) 

where equal and opposite terms are added in equation 
(139). From equation (139), we obtain 


F^Ki[L'~L-(r 2 -l2)]^i(L^L)-Ki(r 2 -~!2) (140) 


where 

L'=ri+lo+l' 2 =histantaneous distance AB 
L'— li+ 10+12= “natural” distance AB with no force 
on the load, i.e. (AB)o 

Using F=F 2 =K 2 (r2-l2), equation (140) becomes 


F = K\(L‘ - L) - Ki 



( 141 ) 


F = +k 2 [u - L] “ Ke * [L “ ^ = K ^ AB - (^)o] 

where K eq is the equivalent stiffness coefficient of the 
springs obtained from 

"ir = ~kr + • 

Equation (141) gives a single expression for the force on 
the load at any instant time. 

The above description presents the best mode con- 
templated in carrying out the invention. The invention 
is, however, susceptible to modifications and alternate 
constructions from the embodiments shown in the 
drawings and described above. Consequently, it is not 
the intention to limit the invention to the particular 
embodiments disclosed. On the contrary, the invention 
is intended to and shall cover all modifications, sizes and 
alternate constructions falling within the spirit and 
scope of the invention, as expressed in the appended 
claims when read in light of the description and draw- 
ings. 

What is claimed is: 

1. A controller for a manipulator wherein the manip- 
ulator and an environment of the manipulator form a 
complex, dynamic model with unknown parameter 
values; the controller comprising: 

an adaptive feedforward controller and an adaptive 
feedback controller including means for applying a 
feedback equation having position and velocity 
control terms; 

said feedforward and feedback controllers being in- 
dependent and multivariable; 

means for generating an auxiliary signal; 

means for summing said auxiliary signal and adaptive 
gains and for outputting a control signal to said 
manipulator; 

said feedforward and feedback controllers having 
adaptive position feedforward and feedback con- 
trol loops, respectively and adaptive position con- 
trol means operating in accordance with the fol- 
lowing position control law: 


Fi=Wi-h) 


( 138 ) 
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O^CO+c^CO+^W+^CO 

wherein R(t) is the desired position trajectory, E^is the 
position tracking error, A, B and C are adaptive feed- 5 
forward gains, and K v are adaptive feedback gains, 
f(t) is an auxiliary signal and F/t) is an applied force. 

2. A controller for a manipulator wherein the manip- 
ulator and an environment of the manipulator form a 
complex, dynamic model with unknown parameter 
values; the controller comprising: 
a feedforward controller and an adaptive feedback ^ 
controller including means for applying a feedback 
equation having force control terms; 
said feedforward and feedback controllers being in- 
dependent and multivariable; 20 

means for generating an auxiliary signal; 
means for summing said auxiliary signal and adaptive 
gains and for outputting a control signal to said 
manipulator; 25 

said feedforward and feedback controllers having 
adaptive contact force control means operating in 
accordance with the following force control law: 


m = PAt) + 



Em + K p {f)E(t) - KAt)Z{t) 


wherein P^(t) is a desired force reference signal, d(t) is 
an auxiliary force signal, KXt), Kp(t) and K v (t) are 
adaptive gains, E(t) is a force tracking error, Z(t) is a 
velocity and F^t) is an applied force. 

3. A controller for a manipulator wherein the manip- 
ulator and an environment of the manipulator form a 
complex, dynamic model with unknown parameter 
values; the controller comprising: 
an adaptive feedforward controller and an adaptive 
feedback controller including means for applying a 
feedback equation having position and velocity 
control terms; 

said feedforward and feedback controllers being in- 
dependent and multivariable; 
means for generating an auxiliary signal; 
means for summing said auxiliary signal and adaptive 
gains and for outputting a control signal to said 
manipulator; 

said controller further comprising position control 
means, force control means; and means for cross- 
coupling said position control means and force 
control means to compensate for said parameter 
values. 
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